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ABSTRACT 

Context. Long-term records of solar radiative output are vital for understanding solar variability and past climate change. 
Measurements of solar irradiance are available for only the last three decades, which calls for reconstructions of this quantity over 
longer time scales using suitable models. 

Aims. We present a physically consistent reconstruction of the total solar irradiance for the Holocene. 

Methods. We extend the SATIRE (Spectral And Total Irradiance REconstruction) models to estimate the evolution of the total (and 
partly spectral) solar irradiance over the Holocene. The basic assumption is that the variations of the solar irradiance are due to the 
evolution of the dark and bright magnetic features on the solar surface. The evolution of the decadally averaged magnetic flux is 
computed from decadal values of cosmogenic isotope concentrations recorded in natural archives employing a series of physics-based 
models connecting the processes from the modulation of the cosmic ray flux in the heliosphere to their record in natural archives. We 
then compute the total solar irradiance (TSI) as a linear combination of the /'' and /'' + 1 decadal values of the open magnetic flux. 
In order to evaluate the uncertainties due to the evolution of the Earth's magnetic dipole moment, we employ four reconstructions of 
the open flux which are based on conceptually different paleomagnetic models. 

Results. Reconstructions of the TSI over the Holocene, each valid for a different paleomagnetic time series, are presented. Our anal- 
ysis suggests that major sources of uncertainty in the TSI in this model are the heritage of the uncertainty of the TSI since 1610 
reconstructed from sunspot data and the uncertainty of the evolution of the Earth's magnetic dipole moment. The analysis of the 
distribution functions of the reconstructed irradiance for the last 3000 years, which is the period that the reconstructions overlap, 
indicates that the estimates based on the virtual axial dipole moment are significantly lower at earlier times than the reconstructions 
based on the virtual dipole moment. We also present a combined reconstruction, which represents our best estimate of total solar 
irradiance for any given time during the Holocene. 

Conclusions. We present the first physics-based reconstruction of the total solar irradiance over the Holocene, which will be of interest 
for studies of climate change over the last 1 1500 years. The reconstruction indicates that the decadally averaged total solar irradiance 
ranges over approximately 1.5 W/m^ from grand maxima to grand minima. 

Key words. Sun: activity; Sun: faculae, plages; Sun: surface magnetism; Sun: solar-terrestrial relations; Sun: sunspots; Sun: UV 
radiation 



1. Introduction 

Radiative energy input from the Sun is the main external source 
of heat to the Earth's coupled atmospheric/oceanic system 
( |Hansen|[2000} |Haigh| 



200 ip. Variations in this external source 



obviously affect the climate, although there is still considerable 
uncertainty and debate on the level of this effect (|H aigh 2001 ; 
[Hansen et al.p002l|2005||Solomon|2007[|Camp & Tung.2007| ). 



One way to disentangle the influence of the Sun on the climate is 
to look for its influence on climate in pre-industrial times, when 
the influence of man-made greenhouse gases can be neglected 
and the natural variations can be isolated. 

Direct measurements of atmospheric and oceanic parame- 
ters during the last century suggest the occurrence of climate 
shifts on several time scales ("Trenberth 1990', Zhang & Song 



2006 Solomon„2007j . Holocene climate proxies also suggest 
the occurrence of millennial scale climate variations ( |Rodbell| 
et al. 1999[ Bond et al. 2001| l. For example, the intensity and 
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frequency of the El Nino and Southern Oscillation (ENSO), 
which is the most important oscillation mode of the atmo- 
sphere/ocean coupled system, seems to have undergone signif- 
icant climatic shifts in the Holocene (Rod bell et aH 1999), with 
strong ENSO events becoming established only about 5,000 
years ago. Furthermore, Holocene climate proxies also suggest 
substantial correlations between atmospheric and oceanic condi- 
tions and solar variations, but the physical mechanisms are not 



understood yet ( [Bond et al.]|2001, Haigh 2001 2007{ [Usoskin 
let al.|2005| ). 

With the current generation of general circulation models it 
is now possible to carry out long-term numerical simulations 
of the evolution of the Earth's climate system (Jungclaus et al. 



2010 Schmidt 2010). Such studies require long time series of 



solar irradiance, ideally covering the whole Holocene. Although 
reconstructions of irradiance since the Maunder Minimum are 



common (e.g. Lean et al. 11995 Solanki & Fligge|l999 Fligge & 



Solanki 2000; Wang et ai.'2005'' Balmaceda '2007';'K rivova et al 
2007 , 2010a; Tapping et al. 2007j]Crouchet al. 2008| l, on longer 
timescales few quantitative reconstructions are available due to 
the lack of sufficiently time-resolved solar data and appropriate 
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techniques for reconstructing the irradiance when only data av- 
eraged over 10 years (i.e. roughly a solar cycle) are available (but 
Steinhilber et al.|2009 1. Here we compute the total (and to a 



see 



limited extend the spectral) irradiance for the Holocene employ- 
ing the reconstructions of the open flux and sunspot number ob- 
tained from the cosmogenic isotope '^C. The motivation for this 



work came from a study carried out by |Balmaceda| ( |2()()7 ), who 
compared the annual and cycle averages of the TSI, total and 
open magnetic fluxes, as computed using the model described 
by [Krivova et al. ( 2007} , and the observed sunspot number. She 
found robust relationships between these quantities. IVieira et al. 
( |200 8) employed these relationships to make a preliminary com- 
putation of the total solar irradiance for the Holocene from the 
sunspot number and open magnetic flux estimated from '^C. 
More recently, Steinhilber et al. ( 2009 1 presented a reconstruc- 
tion of the total solar irradiance for the Holocene based on the 
cosmogenic radionuclide '"Be measured in ice cores. A shorter 
reconstruction has been presented by Delay gue & Bard ( 2010[ ). 

In the present work we derive the relationship between solar 
cycle averaged open magnetic flux and total solar irradiance for 
the physics-based model of [Krivova et al.|(|2007||20 10a| . This is 
then used to reconstruct the irradiance throughout the Holocene 
based on ''^C data. The reconstructed values lying in the tele- 
scope era are compared with the cycle averaged values obtained 
by Krivova et al. ( [2010a ). The present model is so far unique in 
that it is based on a consistently derived, physics-based relation- 
ship between TSI and the magnetic flux, which incidentally does 
not correspond to a direct, instantaneous relationship as is often 
incorrectly assumed. By following this approach the real phys- 
ical processes are considered in all steps, but some parameters, 
which cannot be obtained directly from theory, are evaluated by 
fitting the model's output to observational data. 

We structure the paper as follows. In Sect. 2, we describe the 
physics-based model used to reconstruct the TSI and SSI, which 
is, in principle, identical to the models employed by Krivova 
[et al.| ( [2007[ [2010 a ) to calculate the solar irradiance on time 
scales from days to centuries, but adapted to work with heav- 
ily time-averaged data. In Sect. 3, the reconstruction of the TSI 
and SSI for the Holocene is presented and sources of uncertainty 
are discussed. Finally, we present the conclusions in Sect. 4. 



2. Model Description 

The primary source of information on solar activity in pre- 
telescopic times are concentrations of cosmogenic isotopes, such 
as '^C or '"Be, in natural archives. From these it is possible to 
deduce the Sun's open magnetic flux and the sunspot number, as 
has been demonstrated in numerous publications (e.g., Usoskin 
[erai..2003,.2006aj [Solanki et al. 2004, Steinhilber et al.|2010| l 



In fact, when computing the sunspot number from the open flux, 
the magnetic flux in active regions, ephemeral regions and the 
total magnetic flux are computed. These are also the ingredients 
needed to estimate the total and spectral solar irradiance em- 
ploying the model described by Krivova et al. (2007 2010a[ l. 
The main hurdle is that in general only the cycle-averaged open 
flux is known, so that the models originally developed for time- 
resolved data need to be modified in order to be applicable 
on these long time scales. For the magnetic flux the relevant 
relationships were most recently derived by |Vieira & Solanki 
( 2010| l. Here, we now take the next step and determine the Sun's 
cycle averaged irradiance. 

The aim of this section is to derive the relationship between 
the cycle-averaged TSI and a solar quantity as closely related to 
the observations (isotope concentrations) as possible. The open 



magnetic flux of the Sun is such a quantity. We stress that the 
open magnetic flux itself produces only a rather minor part of 
solar irradiance variations, but is closely linked to the quantities 
which produce larger fractions of the irradiance variations. In 
particular, the cycle averaged open magnetic flux closely follows 
the Sun's total magnetic flux. 

Before we begin, a brief note on terminology: SATIRE 
stands for Spectral and Total Irradiance Reconstructions. The 
basic philosophy of this model family is that all solar irradiance 
variations on time-scales longer than a day are due to the evolu- 
tion of the solar surface magnetic field. Members of the SATIRE 
family of models are SATIRE-S, which models the irradiance 
during the satellite era with high precision (in its present version 
since 1974), and SATIRE-T, which computes irradiance during 
the telescope era (since 1610). Here we introduce the SATIRE- 
M model, where 'M' stands for millennial time scales. 

2. 1 . Model of solar irradiance variations 

The Sun's irradiance responds directly to the amount and dis- 
tribution of magnetic flux on the solar disk, which gives rise 
to dark features, sunspots, and bright features, faculae and net- 
work. In order to quantitatively reproduce solar radiative flux it 
is sufficient to segment the solar disk, using magnetograms and 
brightness images, into dark (sunspot umbra and penumbra) and 
bright (faculae and network) magnetic components as well as the 
unaffected quiet Sun ( Krivova et al.|2003] Wenzler et al.|2007 



2006 j ). The computed irradiance depends both on the area cov- 
ered by a given type of feature and the spatial location on the 
solar disk. In the period before magnetograms became available, 
a simplified version of the model had to be employed. In this 
model a homogeneous distribution of magnetic features over the 
solar surface is assumed, so that only the fractional area of the 
solar disk covered by each of them (global, or disk-averaged fill- 
ing factor, a) and the specific flux, F, (i.e. the flux that would 
be emitted if the whole solar disk was covered by that type of 
feature) of each feature is required. The evolution of radiative 
flux at a given wavelengths. A, can then be computed as follows 
( [Krivova et aT.,2007„2010a| ), 



F(A, t) = aq{t)Fq{A) + a,(t)FM) + ap{t)Fp{A) ■ 
af(t)Ff(A) + a„{t)F„{A) , 



(1) 



where the Fq, F„, Fp, F j, and F„ are the time-independent radia- 
tive fluxes of the five component model: quiet Sun (q), sunspot 
umbra (u), sunspot penumbra (p), faculae (f), and network (n). 
Following Krivova et al. ( 2007[ ) and others, the spectrum of each 
component that was calculated by [Unruh et al.| ( 1999 ) using the 
ATLAS9 code of Kumcz (1993). Here, we employ the same 
model to describe the brightness of facular and network com- 
ponents. 

The fraction of the solar disk covered by each component is 
represented by Oq, a„, a^, aj and a„, respectively. By taking this 
approach we are assuming a nearly homogeneous distribution of 
magnetic features on the solar disk. The filling factor of the quiet 
Sun can then be written as: 



aq{t) = 1 - au(t) - Upit) - a fit) - a„(t) . 
Substituting Eq. (j2|i in ([T]), we obtain 



(2) 



F(A, = (1 - a,.(t) - ap(t) - af(t) - a„(t))Fq(A) + a,(t)F,XA) + 
apit)Fp(A) + afFfiA) + a„(t)F„(A) . (3) 
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This equation can be rewritten as 



and the constant c is given by 



F(A,t) = F^(A) + aMAF„(A) + apmFp{A) + 

afmFf(A) + a„(t)AF„{A) , (4) 



where 

^Fu^pj^„(A) = Fi,^pj\„(A) - Fq(A) . 



(5) 



2.2. Filling factors for sunspot umbrae and penumbrae 

We compute the filhng factors for the umbral and penumbral 



components by employing the approach described by Krivova 
et al. (2007V A fixed ratio of 0.2 between umbral and total spot 



area is assumed following |Wenzler et al., ( ,2005| ), which translates 
into 



Qu/a^ = aj(a„ + a„) = 0.2 . 



(6) 



We compute the sunspot area, i.e. the fraction of the disk cov 
ered by all sunspots on the solar disk, by making use of a linear 
relationship to the sunspot number {R) ( FUgge & Solanki|1997 



Balmaceda et al.|2009[[Hathaway|2010 



Us - AiR + A2 , 



(7) 



where Ai is the proportionality coefficient and A2 is the offset. 
The key step for extending the irradiance model to the whole 
Holocene is to compute the sunspot number as well as the mag- 
netic flux in active regions (AR) and ephemeral regions (ER) 
from the open magnetic flux, since it is the first solar parame- 
ter obtained from the chain of models that convert cosmogenic 
isotopes data into solar parameters. This problem was first ad- 
dressed by Usoskin et al. ("2002 1 and later updated by Usoskin' 
|eral.,(2003 2007) , Solanki eTaL] ( |2004| l, and Vieira & Solanki, 
(|2010|l. liTttie cur rent work we employ the method described by 
Vieira & Solanki (2010') to compute the total solar surface mag- 
netic flux and the sunspot number from the open flux. From the 
information recovered from cos mogenic isoto pes, in particular 
the INTCAL98 ''^C record (Stuiveretal.|l998l, we are limited to 



estimate 10-year averages of the sunspot area and the solar sur- 
face magnetic components. Following Vieira & Solanki ( 2010| l 
we write for decadal values 



< R >j- flR < 



'open 



>j +bR < 



^open '^y+l 



(8) 



where < ... >j and < ... >y+i are the / and / + 1 10-year 
averaged values (or data points) of a variable. Since 10-years is 
close to the average length of the solar cycle (approximately 1 1 
years for the past 23 cycles), and + 1 also roughly refer to 
the and /'' + 1 cycles. The coefficients qr and bR are given by 



Qr 



and 



1 

CTl 



cAf 



(9) 



(10) 



where Ar is the samphng interval, t\ is given by 



open 



TephkX 



eph 



max,2 1 
l>max.2\ 



(12) 



In Eqs. ( (8p2| ) 



Tgp^„ and Tgp^i^ are the decay times for the rapidly 



eph 



and slowfy evolving components of the open flux. t;J^, and r 
are the flux transfer times from active and ephemeral regions to 
the slowly evolving open magnetic flux, while tJJ^, is the flux 
transfer times from active regions to the rapidly evolving open 
magnetic flux, e™^'^' and R'™"'^'^ are the maximum emergence 
rate of magnetic flux in active regions and the maximum value 
of the sunspot number observed in cycle 21, respectively. The 
parameter k is the ratio between sunspot number amplitude and 
the 10-year averaged value ( [Usoskin et al.|2007 1. X is the scal- 
ing factor for the flux emergence rate in ER relative to active 
regions. Here, we employ the set of parameter values employed 
by Krivo va et al.| ( 2010a) (see Table 1). The sampling interval, 
Af, is 10 years. The effective decay time scales of active regions 
(AR), Tact, and ephemeral regions (ER), Tgph, are given by 



1 

Tact 

and 
1 



1 1 1 

-77- + —, 1 



act act 



1 1 

77— + 



(13) 



(14) 



'^'■P'' '^eph '^eph 

where, r^^, and t°^,j are the decay time scales of active and 
ephemeral regions, respectively. 

Combining Eqs. (|6]), (|7]) and (|8]), we obtain 



< a„ > i 



and 



0.2(AiflR < 



^open 

>j +AibR < 



+A2) 



a,, < 



Popen -^J -^bii < I 



(15) 



< ap >j = O.^iAiOR < cf)open >j +AibR < (popen > j+\ +A2) 

1^0 pen 

>j+\+Cp. (16) 



^open -^bp < 



The sunspots of a given cycle influence not only the open 
flux of that particular cycle, but also of the next cycle due to the 
slow decay of the open flux ( Solanki et al.|2000 ). Consequently, 
it takes the knowledge of the open flux from 2 cycles to get the 
sunspot areas and the magnetic fluxes of active and ephemeral 
regions for one cycle, as indicated by Eqs (15i and (16). It 



means that the part of the magnetic flux that emerges in ac- 
tive and ephemeral regions, and is dragged outward by the solar 
plasma and permeates the heliosphere, is implicitly represented 
in Equation (js]). The relative contribution between the and 
+ 1 data points (aR/bR = (Af - T^pe„)/-r;',^,„,) depends on the 
decay time scale of the slow component of the open flux and on 
the sampling interval. For the set of parameters employed here 
(given in Table 1), the contribution of the data point relative 
to that of the -i- 1 data point is 



Or 
bR 



b„ 



2.4, 



(17) 



1 



1 



1 

Af 



(11) 



which suggests that a direct, instantaneous relation between 
sunspot areas and the open flux is not appropriate. 
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2.3. Filling factors for faculae and networi<. 

We compute the 10-year averaged filling factors for faculae 
(< af >j) and the network (< a„ >j) from the solar surface 
magnetic flux employing the same scheme as in previous appli- 
cations of the SATIRE mod els (|Krivova et al.|2003H2007||2010a] 
Wenzler et al.|2004j,2005, , 2006| . <af >j and < g„ >j are pro- 
portional to the magnetic flux m faculae (< cpf >j) and network 
(< (p„ >j), respectively, until a saturation limit (ps„, is reached. 
The saturation flux is a free parameter of the SATIRE models, 
which depends on the resolution and the noise level of the em- 
ployed magnetic field maps. Observing that the 10-year averaged 
values < <pf >j and < (p„ >j are always smaller than the satura- 
tion level, so that they never saturate, we can write 



The contribution from the next i.e. (j + l)-st cycle is due to the 
extended period over which ephemeral regions emerge. From 
Eqs. (|24| and ( |25| ), 

(26) 



< 



where 



^eph ' 



-kXoR + 1 



and 



(,21 



Teph - 



kXbR 



' J^max,21 

Substituting Eq. (|26 



into 



19 1, we obtain 



(27) 



(28) 



1 



and 



< a„ > i 



sun ^ sal, f 



sun ^ sal, n 



(18) 



where S „,„ is the Sun's surface area. Bsaij and Bsat.n are the satu- 
ration field strengths for faculae and network, respectively. Here 
we follow [Krivova et aT] ( |2007| and [Krivova et al.|P010a[ ) who 
incorporated a correction in Bsaij that takes into account that the 
faculae mainly appear in the active belts. 

The magnetic flux in active regions (< (pac >j) is the sum of 
the flux in sunspots (umbrae plus penumbrae) and faculae. Thus, 
the magnetic flux in faculae is given by 

< >, = < (pac, >i - < <Pu >i - <'l>p >i ■ (20) 



t'open 



< 



fl„ < 



'open 



S simBsat,n 

>j +b„ < 



(29) 



(19) While the relative contributions between and + 1 data 
points of the open flux for the umbral, penumbral, and facular 
filling factors depend on the decay time of slowly evolving open 
flux and the sampling time, the relative contribution to the net- 
work is given by 



b„ 



Following Krivova et al. ( 2007| l, we compute 

< <Pu >j + <<Pp >r< au >i [Bz,u\Ssun+ < ap >j [Bz,p]S sun, (21) 

where [B^.,,] and {B~p] are the area-averaged vertical compo- 
nents of the magnetic field strength in umbrae and penumbrae, 
respectively. Note that since < Ou >j and < Op >j are indirectly 
obtained from < R >j, it is assumed implicitely that averaged 
over 10 years there is no major difference between the number 
of sunspots on the visible part of the solar surface (solar disk) 
and on the hidden part. 

From Eq. (25) of |Vieira & Solanki| ( |2010| l, the magnetic flux 
in active regions is approximately given by 



Af 



' open 



eph' 



kX 



1 



(30) 



For the set of parameters employed here, this ratio is approx- 
imately 3.3. This indicates that although the contribution from 
the next cycle is lower than that from the cycle for which we 
want to know < a„ >y, it is still necessary to know the next state 
in order to evaluate the filling factor of the network accurately. 

2.4. Estimate of the total solar irradiance 

The solar irradiance at wavelength A can be computed by sub- 
stituting the expressions for the filling factors obtained for the 
components into Eq. Q 



< F(A) >j= OfiA) < . 
where 



>j +bf{A) < < 



'open ^7+1 



+F,(A), (31) 



apiA) = auAFu(A) + a,,AFp{A) + afAFf(A) + a„AF„(A) , 
and 



(32) 



' J^max,2\ 



(22) 



Combining Eqs. ( 15 i, (18 1, (19 1, (20 1, (21 1, and (E2b, we ob 



tain 



< Of >j- Of < (popen >j +bf < (popen >>+! +Cf . 



(23) 



The sum of ephemeral regions flux and of the open flux de- 
scribes the evolution of the network: 



eph 



(24) 



From Eq. (29) of Vieira & Solanki ( 20I0| l, the lO-yeai^ averaged 
magnetic flux in ephemeral regions is 



beph >j 



"^eph 



^eph 



^-kX<R>i 

max, 2 1 ■' 



R' 

max,2 1 

^^^kxL 

J^max,2\ \ 



Or < ' 



>j +bR < (popen >i 



bpiA) = b,AFu(A) + bpAFpiA) + bfAFf(A) + b„AF„{A) ■ (33) 

In order to obtain the total solar irradiance one can either in- 
tegrate < F{A) >j over all wavelengths or right from the begin- 
ning work with wavelength integrated radiative fluxes and flux 
differences 



Jo 



Fg(A)dA, 



and 



AF' f 



[Fu,pj,n{A)-F,(A)\dA. 



(34) 



(35) 



Replacing the Fq and AFj in Eqs. (pT|i-(33 i by F'^ and AF'., 
i = u, p, f, n, one then obtains the total solar irradiance averaged 
over 10 years, < F' >j. We refer for clarity to < F' >j and 
< F'l >j, i = M, p, f, n, as TSI and TSIj. 

Hence, the solar spectral and total irradiance can be com- 
puted from the and /'' + 1 values of the open flux in a rela- 
tively straightforward manner. 
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2.5. Reconstruction of the open flux from cosmogenic data 

The open flux in the past can be reconstructed from data of cos- 
mogenic isotopes measured in terrestrial archives. Here we use 
data from ''^Cmeasured in tree-rings for the past millennia as 
compiled by the INTCAL04 coUaboration (Rei mer et al.|2004| l. 
We have applied a commonly used approach of inverting the cos- 
mogenic isotope data into cosmic ray flux. This approach, pro- 



posed Usoskin et al. (2003 2004j|, is now a standard method 
used to reconstruct of the solar or heliospheric parameters in 
the past (e.g. "Solanki et al. 2Q0 4| [Usoskin et al. 2006a, 2007; 
Vonmoos et al. 1,2006, iMuscheler et al.||2007[ [Steinhilber et al.| 
2008 1. The method includes the following steps. 



1. First, the measured values of A''*C are converted, using a 
carbon cycle model, into the global '"*C production rate, 
Qi4c. We used a standard multi-box carbon cycle model 
([Siegenthaler 1980| l parameterized using the approach of 
[Usoskin & Kromer| ( [^05l ). 

2. Considering prescribed geomagnetic fiel d data and a model 
of '"*C production in the atmosphere (Castagnoli & Lai 
i980j ), the production rate Q''*C was converted into the cos- 
mic ray spectrum quantified via the modulation potential O. 

3. Using a standard theory of cosmic ray transport in the 
Heliosphere one can relate the modulation potential O to the 
open magnetic flux assuming a constant 1 1-yr averaged solar 
wind speed (e.g. [Usoskin et al.|2002j ). 

Here we give a brief description of the method, while full details 
are provided elsewhere ( Solanki et al. 2004 Usoskin et al. 2007| l. 

The uncertainties of the reconstructed open flux cr^^^^^ have 
several sources. The major source of the systematic error is re- 
lated to the uncertainties of the '"^C production rate in the pre- 
Holocene period. Other sources of systematic uncertainties are 
related to the paleomagnetic models (considered below in de- 
tail), uncertainties of the carbon cycle model, uncertainties of 
the ^'^C production model, and uncertainties of the cosmic ray 
heliospheric transport model. Random errors are related to mea- 
surement errors in A''*C and in geomagnetic data. All these 
uncertainties are described in detail by S olanki et al.[ ( [2004| 
Supplementary matarial). The total uncertainty in the open flux 
•^"Aop™ is °f '^he order of 0.2 ■ 10'"* Wb for the last five millennia 
with relatively good geomagnetic data and rises up to 0.5 ■ 10'"* 
Wb in the early Holocene. 

We note that computation of the open flux from '^C is not 
a part of the present work as it was performed using the same 



algorithms and codes as in Usoskin et al. ( 2007 1, but now ap- 
plying different paleomagnetic models. We point out that the 
14C-based open flux is consistent with the '^'^Ti data from mete- 
orites ( [Usoskin et al.|2006bl ). 

3. Results 

3.1. Model setup and validation 

The equations modeling the evolution of TSI during the 
Holocene have been derived in the previous section. In this sec- 
tion, we present some tests. We compute the coefficients for the 
solar spectral and total irradiance model for the Holocene ac- 
cording to the parameters of the SATIRE-T model presented by 
Krivova et al. ( 2010a[ l. Table|2]lists the coefficients entering the 
equations given in Sect. 2 that are used to compute the solar irra- 
diance from the open magnetic flux (which is deduced from '"^C 
data). We stress that we have not introduced new free parame- 
ters when simplifying the model to work with 10-year averaged 



values. Thus, uncertainties in the parameters are a heritage from 
the SATIRE-T model described by [Krivova eTaLl ( |2010a[ ). 

[Krivova et al. (20I0ai constrained the parameters of the 
model by comparing the model output with observations and 
reconstructions of five solar variables. While two data sets are 
magnetic flux variables, three data sets are related directly to the 
solar irradiance. The total solar surface magnetic flux was con- 
strained by comparing the model output with estimates carried 
out at the Mt. Wilson Solar Observatory (MWO), National Solar 
Observatory Kitt Peak (KP NSO) and W ilcox S olar Observatoiy 
(WSO) over cycles 20 to 23 ^Arge et al. 200 2; Wang et al.[2005| l. 
The open flux was compared to the reconstruction based on the 
geomagnetic aa-index ( jLockwood et al.|2009a|b[ ) from 1904 to 
2008. These two data sets provided constraints on the ampli- 
tude of the magnetic activity cycle and its long-term evolution. 
We stress that the evolution of the open flux as reconstructed 
based on the aa-index, which nearly doubled during the last 
century, is well reproduced by the magnetic field model (see 
Figure [T^). In addition to the reconstruction based on the geo- 
magnetic aa-index, the long-term evolution of solar magnetic ac- 
tivity can also be constrained by cosmogenic isotopes deposited 
in terrestrial and extra-terrestrial archives (Sol anki et al.[[20()4| 



Usoskin et al. 2006b i. While terrestrial archives can span the 



entire Holocene, the processes leading to the deposition of the 
cosmogenic isotopes are affected by the evolution of, e.g., the 
Earth's magnetic field and of the highly coupled atmosphere- 
ocean system. Extraterrestrial archives, such as meteorites, are 
not affected by such processes and can be employed to test the 
long-term trend observed in the modeled open flux. In order to 
substantiate the analysis carried out by Krivova et al. ( 2010a[ l 
and Vieira & Solanki (2010) , we test the reconstruction of the 
open flux by comparing the activity of cosmogenic isotope "''^Ti 
in meteorites, which fell during the past 240 years, and the mod- 
eled activity based on the open flux computed in the framework 
of the above models. We carried out the analysis using the same 
method as described by Usoskin et al. (2006b) including evalu- 
ation of all known sources of uncertainties. Figure [T]5 presents 
a comparison between the ''^Ti measurements in 19 meteorites, 
indicated by diamonds, and the modeled '^"'Ti activity at the time 
of the meteorite's fal l (see details in Uso skin et al.| ( |2006b| l and 



Taricco et al. 



( [2006| l). The '♦^Ti activity is given in disintegra- 
tions per minute per kilogram (dpm kg''). The is about 7.6 
for the 19 data points, i.e. about 0.4 per degree of freedom, that is 
comparable to the results previously discussed by Usoskin et al. 
( [2006b I. Thus, the open flux modeled employing the parameters 
of the SATIRE-T model follows the data reasonably well, which 
leads us to conclude that the parameterization used is consistent 
with the available longer term solar data. 

The model discussed by [Krivova et al.| ( [20I0a) l was fur- 
ther constrained by requiring the computed TSI variations to 
match the PMOD composite (IFrohhch 20061, which is ob- 
tained from space -based observations since 1978. Here, we em- 
ploy (among other measurements) the difference between the 
SATIRE-T model output and the PMOD composite to estimate 
the uncertainty of the TSI model for the Holocene. We refer to 
this difference as the model error. Figure |2^ presents a com- 
parison between the observed (PMOD composite; blue line) 
and the modeled (SATIRE-T; green line) TSI. The plot cov- 
ers the period from Nov./1978 to Dec./2003, which was em- 
ployed by Krivova et al. ( 2010a| l to constrain the SATIRE-T 
model. Note that as the absolute value of the TSI is not well 
determined we offset the model to match the PMOD compos- 
ite level, although observations employing the SORCE/TIM in- 
strument suggest that the absolute value during the minimum 



5 



Vieira, Solanki, Krivova, and Usoskin: Evolution of the TSI/SSI 



between cycles 23-24 is approximately 4-5 W/m^ lower than 
that measured by other current instruments (Kopp et al.||2005 



Kopp & Lean|20H ). Figure |2|3 presents the difference between 



the PMOD composite and the adjusted SATIRE-T model. The 
gray and black lines display daily values and 90-day means, re- 
spectively. Figure|2j; presents the distribution of the model eiTor 
The standard deviation (cr) and the median absolute deviation 
(MAD) are indicated in the panel. The MAD value is defined 
as median(abs{X - median{X))), where X is the sample popu- 
lation. The red line represents the corresponding normal distri- 
bution for the estimated mean and standard deviation. While the 
empirical cumulative distribution function (CDF) of the model 
error is presented in Figure [2}l, the probability plot for a Normal 
Distribution is presented in Figure|2^. In the probability plot, the 
marks ('x') indicate the empirical probability versus the model 
error values. Note that the scale of the ordinate axis is not lin- 
ear The distance between the ordinate tick marks is proportional 
to the distance between the quantiles of a normal distribution. 
The dashed line connects the 25''' and the 75''' quantiles. The as- 
sumption of normality is reasonable if the data points are near 
the dashed line. However, this plot clearly indicates that the dis- 
tribution is not normal. This is a consequence of the procedure 
employed to constrain the model, which reduced the weight of 
outliers to compute the x^- Here, we use the MAD, which is 
a more robust estimator than the standard deviation and vari- 
ance, as a measure of the dispersion of the model error distri- 
bution. In order to estimate the dispersion from the MAD value 
we multiply the MAD value by the scale parameter for the nor- 
mal distribution ( Mosteller & Tukey|1977 jl. The value obtained 
is approximately 0.28 Wlnt^. Note that this uncertainty applies 
to daily data, whereas we are dealing with decadal averages. For 
these the difference between reconstruction and observations is 
very much smaller. Nonetheless, we retain 0.28 W/m^ as the 
model error, since it corresponds roughly to the uncertainty in 
the Maunder Minimum level (Krivova et al. 2007). Figure |2]^ 
presents the scatter plot of the observed (PMOD Composite) ver- 
sus modeled TSI (SATIRE-T). The green line presents the best 
linear fit to the data set. The coefficients and their uncertainties 
are indicated in the figure as well as the correlation coefficient 
(Rc = 0.79). 

We validate the SATIRE-M model during the telescope era 
by comparing its output with that of the SATIRE-T model. The 
two models are not independent, so that with these comparisons 
we mainly test the simplifications introduced in order to be able 
to adequately deal with the 10-year averaged data. To compare 
both models, first we compute 10-year averages from 1700 AD 
to 2000 AD of the open flux, total magnetic flux and TSI as 
well as of the contributions from faculae, network and sunspots 
(umbrae -i- penumbrae) to TSI from the SATIRE-T model. We 
then employ the 10-year averaged open flux from SATIRE-T as 
input to SATIRE-M. This allows us to compute the contribution 
of each atmospheric component to TSI employing SATIRE-M. 

Figure |3] presents a comparison between the output of the 
models SATIRE-T and SATIRE-M. Each panel contains the fol- 
lowing elements: The red lines indicate the set of expectation 
values for the model, while the blue lines are the best fit lines. 
The linear model coefficients and the correlation coefficient are 
also shown in each panel. The vertical error bars correspond to 
the dispersion of the SATIRE-T model error computed from the 
MAD value of the SATIRE-T error distribution and corrected 
for the normal distribution. The estimate of the SATIRE-M er- 
ror is discussed in the next subsection. In addition to the disper- 
sion of the TSI error distribution, the other constraints imposed 
by |Krivova et al.| ( |2010a| l can also be employed to compute the 



dispersion of the model error distributions of the total magnetic 
flux and components of the TSI. Additionally to the constraint 
imposed by the PMOD composite, [Krivova et al. (2010a) con- 
strained the amplitude of the facular contribution to TSI by com- 



paring the model output with the estimate computed by Wenzler 



et al. (2006). The dispersion of the difference between the mod- 
els is about 0.27 W/m^ when corrected for the normal distribu- 



tion. Although the sunspot deficit computed by Wenzler et al. 



( |2006 1 was not employed to constrain the SATIRE-T model, it 
could be employed to estimate the uncertainty of the sunspot 
(umbral + penumbral) contribution to TSI. In this case, the dis- 
persion is about 0.14 W/m^. In order to ensure that not only the 
total UTadiance is reproduced correctly, we also compare the so- 
lar irradiance integrated over the wavelength range 220-240 nm 
(see'Krivova et al.'2009V2010bi) with the SATIRE-T model out- 
put by Krivova et al. (2010a). 

The SATIRE-M model reproduces quite accurately the evo- 
lution of the total magnetic flux as computed by the SATIRE- 
T model {Rc = 0.98; slope = 1.00 + 0.02; Figure [3]), but 
it slightly overestimates the magnitude of the total magnetic 
flux {offset = -1.88 ± 0.79 (xlO"* Wb)). The SATIRE-M 
model also reproduces well the total solar irradiance computed 
employing SATIRE-T {Rc - 0.98), but slightly overestimates 
low values and underestimates high values, which result in a 
slope - 1.09 ± 0.05). A similar behavior is also observed when 
comparing the facular and network brightness excess. The co- 
efficient of correlation is approximately 0.98 while the slope is 
approximately 1.05. The sunspot (umbrae + penumbrae) bright- 
ness deficit is also well reproduced but the coefficient of corre- 
lation {Rc - 0.96) is lower than the one obtained for the com- 
parison of the brightness excess from faculae and network. The 
slope of the linear regression is close to unity {slope - 0.98). 
The UV irradiance integrated from 220 to 240 nm follows the 
behavior of the facular and network contribution of the TSI, with 
a high coefficient of correlation {Rc = 0.99). 

The correlation coefficients for all cases exceed 0.96 and the 
differences between the models display nearly Gaussian distri- 
butions. Thus, we conclude that the 10-year averages of the TSI 
(and UV iiTadiance) can be obtained from a linear combination 
of the /'' and /'' + 1 observations of the open magnetic flux with 
reasonable accuracy. 

3.2. Reconstructions of the solar irradiance 

As shown in Sect. 2, an estimate of the evolution of the TSI 
during the Holocene can be obtained from a linear combination 
of the /'' and /'' + 1 decadal values of the open flux, which 
can be computed from cosmogenic isotopes recorded in natu- 
ral archives. As described in Sect. |2.5| the reconstructions of 
the open magn etic flux used here ar e based on the INTCAL04 
A C data set (Reimer et al. 112004') employing almost exactly 
the same algorithm as the one used by Solan ki et al.| ( |2004| 
and [Usoskin et al.] ( |2007| [2009[ ). In addition to the uncertain- 
ties in computing the TSI from decadal values of the open flux 
employing SATIRE-T, there are also uncertainties related to the 
computation of the open flux from cosmogenic isotopes (see 
Sect. 2.5 I. The largest uncertainty is posed by the evolution of 
the geomagnetic field ( [Usoskin et al.|2006a 2007, Sno wball &[ 
Muscheler,2007) ) . Information of the geomagnetic field prior to 
the era of direct measurements can be retrieved from archeo- 
logical artifacts, lava flows, and sediments (see [Donadini et al.[ 
2010, and references therein). In order to reconstruct the evo- 
lution of the geomagnetic field it is necessary to determine at 
least one component of the magnetic field and the period that the 
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Fig. 1. (a) Evolution of the solar open magnetic flux. The black and red lines present the estimates of the open flux from the 
geomagentic aa-index and from observations of the interplanetary magnetic field ( [Lockwood et al.|2009a|b )l, respectively. The blue 



line display s the modeled open flux from SATIRE-T. (b) Ti activity in stony meteorites as a function of the time of fall (see Taricco 
|et al .'200"6^ and references therein). Er ror bars refer to 1-sigm a uncertainties. The red line corresponds to the modeled ^"*Ti based on 
the SATIRE-T open flux according to Usoskin et al. (2006b I. The '^'^Ti activity is given in disintegrations per minute per kilogram 
(dpm kg ') of Ni + Fe. The dashed red hnes correspond to the estimated 1-sigma uncertanty of the model 



information was stored in the natural archive. Both parameters 
introduce uncertainties to the determination of the Earth's mag- 
netic field, as discussed bylDonadini et al. (20101). Paleomagnetic 



reconstructions of the Earth's dipole moment can be based on 
different assumptions. Precise determination of the true dipole 
moment is hardly possible before the era of systematic geomag- 
netic measurements started in the 19th century. When the quality 
and spread of data is still good enough, a VDM (virtual dipole 
moment) can be estimated, assuming that the geomagnetic field 
can be approximated by a tilted geocentric dipole. When the 
spatial coverage of paleomagnetic samples becomes too patchy, 
only the VADM (virtual axial dipole moment) can be estimated 
by assuming that the geomagnetic axis is aligned with the geo- 
graphical axis. Both VDM and VADM form a proxy for the true 
dipole moment (DM), but they include also non-dipole contribu- 
tions assigned to the dipole strength. Since the VADM assumes 
the dipole to be aligned with the geographical axis, its values 
can be strongly distorted if the dipole axis is tilted. VADM and 
VDM estimates may differ significantly from the true DM, es- 
pecially on short time scales due to the non-dipole influence in 
the VADMs caused by a strong regional bias in data distribution 
( [Korte & Constable|20"05] l. Generally VDM models give a better 
approximation to the true DM when the data coverage is suffi- 



cient but tend to underestimate it with spatially limited data sets. 
In contrast, VADM models tend to overestimate the true mo- 
ment, especially for insufficient data coverage and tilted dipole 
axis ( Donadini et al.|[20 10). Therefore, we apply the following 
approach to estimate the systematic uncertainty of the results due 
to paleomagnetic reconstructions: we consider the VDM to be a 
lower bound and VADM an upper bound to the true dipole mo- 
ment. In order to be able to provide a single reconstruction of TSI 
over a certain period of time, we simply average the output from 
VDM and VADM based reconstructions of TSI. We do not claim 
that the average of the two groups is the most probable value, al- 
though the most probable value lies between these bounds. Table 
3 presents the geomagnetic models employed here, references, 
and the time coverage. Note that we have employed the global 
reconstructions from penevey et al.| ( |2008) . Figure |4] displays 
the reconstructions of the geomagnetic field intensity, including 
VADM reconstructions by |Knudsen et al.| ( |2008|l and|Genevey 
et al. ( 2008[l, and the VDM rec onstructions by Genevey et al. 



( |2008| l and [korte & Constable] ( [20051 )). Here, we consider this 
set of reconstructions to estimate the uncertainty of the TSI re- 
construction caused by uncertainties in the geomagnetic field. 
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Fig. 2. (a) Comparison between the observed (PMOD Composite; blue line) and the modeled (SATIRE-T; green line) TSI from 
Nov/1978 to Dec/2003, (b) Difference between the observed and modeled TSI. The gray and black lines present the daily and 
90-day averaged values, respectively, (c) Distribution of the difference between daily values of the observed and modeled TSI 
(model error). The standard deviation and the median absolute deviation are indicated in the panel. The red curve represents the 
corresponding probability distribution function for the normal distribution with the mean and standard deviation of the model error, 
(d) Cumulative distribution function (CDF) for the difference between observed and modeled TSI. (e) Probability plot for normal 
distribution with the mean and standard deviation of the model error, (f) Scatter plot of the observed versus modeled TSI. The red 
line indicates the set of expectation values of the model. The correlation coefficient and the coefficients of the best linear fitting 
(green line) are presented in the panel. The p-value for testing the hypothesis of no correlation is lower than 0.5, which indicates 
that the con-elation is significant. 



3.2.1 . Reconstructions for the last 3,000 years 

In Figure|5] we present the reconstructions of the TSI for the last 
3,000 years, which is the period over which all reconstructions 
of the geomagnetic field overlap. Panel (a) presents the recon- 



structions of the TSI. The reconstructions based on the VADM 
models by Knudsen et al. (2008 1 and Genevey et al. (2008| are 
presented as blue and thin red lines, respectively. The reconstruc- 
tions based on the VDM models by Korte & Constable] ( |2005| l 
and |Genevey et al. (2008jl are displayed as the gray and cyan 
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y = (0.98 ± 0.06) * x + (0.03 ± 0.02) 
R = 0.96 (p-value < 0.5) 
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Satire-M UV (10° VJ/m" nm) 



Fig. 3. Comparison between 10-year averaged values of the parameters computed employing the SATIRE-M and SATIRE-T models 
from 1700 AD to 2000 AD. The following scatter plots are presented; (a) Total magnetic flux; (b) Total solar irradiance; (c) con- 
tribution of the bright components (faculae -i- network) to the total solar irradiance; (d) Sunspot (umbra + penumbra) contribution 
to TSI; (e) UV integrated in the range 220-240 nm. The red lines present the set of expectation values for the models and the blue 
lines are the best linear fits. The correlation coefficient and the coefficients of the linear regression are presented in the panels. 



lines, respectively. This color scheme applies for panels (a), (b) 
and (d). For reference, the SATIRE-T reconstruction is presented 
from 1700 to present (thick red line). The different VADM mod- 
els give results very similar to each other over the full plotted 
period of time, whereas the VDM models give somewhat higher 



irradiances, in particular the VDM reconstruction based on Korte 



& Constable ( 2005 ), because of their systematically lower pale- 
omagnetic dipole moment (see Fig. |4] Sect. |3.2| and Appendix 
[Al. 



Uncertainties in the reconstructions are presented in Fig.|5]3- 
d. Here, we take into account three major sources of uncer- 
tainties to the TSI computed from 10-year averaged values of 
open flux. The first one is the uncertainty in the TSI computed 
employing the SATIRE-T model {cttsi-satire), which is rep- 
resented by the red dotted lines in panels (b) and (c) and was 
discussed in Sect. 3.1. The second source of uncertainty lies in 
the computation of the open magnetic flux from cosmogenic iso- 
topes (o"0,^„,), which propagates into the TSI employing the pro- 
cedure described in Sect. 2. Panel (b) presents the estimate of the 
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Table 1. Parameters of the TSI model presented by Krivova et al. (2010a i and employed here. 



Parameter 


Symbol 


Value 


Unit 


AR Flux decay time scale 




0.30 


year 


AR Flux to Slow Open Flux transfer time scale 




71.2 


year 


AR Flux to Rapid Open Flux transfer time scale 


' act 


2.1 


year 


ER Flux decay time scale 




eph 


0.0016 


year 


ER Flux to Slow Open Flux transfer time scale 


t" , 

epn 


17.8 


year 


Rapid Open Flux decay time scale 


t'' 

open 


0.16 


year 


Slow Open Flux decay time scale 


open 


2.97 


year 


ER amplitude factor 


X 


78 




ER cycle extension parameter 


Cx 


5.01 




Saturation flux in faculae 


4>satJ 


260.2 


G 


Saturation flux in network 




800 


G 


Umbral brightness deficit 


A/„ 


-871.8 


Wm.-^ 


Penumbral brightness deficit 


Alp 


-277.1 


Wm.-^ 


Facular brightness excess 


Mf 


85.6 


Wm-^ 


Network brightness excess 


AI„ 


85.6 





Table 2. Coefficients of the model to compute TSI for the Holocene (SATIRE-M). The dependent variables (listed under the heading 
'Quantity') are computed from the equation x = a < (popen >j +b < <popen >j+\ +c 



Quantity 


a 




b 




c 






a/b 


<Pac, 


2.4 




1 










2.37 


4>eph 


2.5 




1.1 










2.37 


a,, 


2.79e-005 




1.17e-005 


[Wb-'] 


-9.3e-007 






2.37 




l.lle-004 




4.70e-005 


[Wb-'] 


-3.7e-006 






2.37 


«/ 


1.99e-003 


[wr'] 


8.39e-004 


[Wb-'] 


2.6e-005 






2.37 


a„ 


7.18e-004 


[Wb-^] 


2.16e-004 


[Wb-'] 








3.32 


TSI, 


-0.024 


[Wm-^Wb-^] 


-0.01 


IWm-^Wb-'] 


8.14e-004 


[Wm 




2.37 


TSIp 


-0.031 


[Wm-^Wb-^] 


-0.013 


[Wm-^Wb-^] 


1.04e-003 


[Wm 




2.37 


TSIf 


0.1 


[Wm-^Wb-^] 


0.043 


[Wm-^Wb-^] 


1.37e-003 


[Wm 




2.37 


TSI,, 


0.061 


[Wm-^Wb-^] 


0.019 


[Wm-^Wb-^] 








3.32 


TSI 


0.11 


[Wm-^Wb-^] 


0.039 


IWm-^Wb-'] 


1364.72 


[Wm 




2.83 


(/V[220-240 nm] 


2e+006 


[Wm-^nmWb-'] 


7.4e+005 


[Wm-^nmWh-^] 


9e+008 


[Wm-^ 


nm] 


2.67 



Table 3. Paleogeomagnetic Models 



Model 


Reference 


Period Covered 


Concept 


KN08 


(1) 


9495 BC - 1705 AD 


VADM 


GN08-8k 


(2) 


5995 BC - 1705 AD 


VADM 


GN08 


(2) 


995 BC - 1705 AD 


VDM 


KC05 


(3) 


4995 BC - 1705 AD 


VDM 



References. (1) 


Knudsen et al.|12008|l 


|Korte & Constab 


le|(|2005|l 



uncertainty of the j data point due to the open flux (o"j5/_0 - 

(ars/0-^„p„,)^ + (bTsio-^„,,y)- The third source of uncertainty 
is the discrepancy between the reconstructions based on con- 
ceptually different paleomagnetic models {(Ttsi-cm)- Panel (c) 
shows the difference between the maximum and minimum val- 
ues of the TSI value computed employing different paleo- 
geomagnetic models {ctjsi-gm)- Panel (d) presents the total un- 
certainty that is given by the sum of the evaluated uncertain- 
ties = (rlsi-SATiRE + '^rs/-*,,,™ + (^Isi-cm)- For this pe- 
riod, the major sources of uncertainty are related to the irradi- 
ance model itself and to the conceptually different paleomag- 
netic models. The uncertainty introduced by the open flux is 
much lower Moreover, the overall uncertainty during this period 
does not exceed 0.4 W/m^. Note that we have followed common 
practice and have assumed that there were no significant changes 



in ocean mixing during the Holocene ( Stuiver|1991 Usoskin & 
|Kromer|20 05l. 

A summary of the statistical properties of the TSI reconstruc- 
tions for the last 3000 years is plotted in Fig. A. 1 of the appendix. 



From this figure and statistical tests we found that the distribu- 
tions do not correspond to normal distributions. Additionally, re- 
constructions from the same group (VADM or VDM) have the 
same median values. However, the reconstructions belonging to 
different groups disagree from each other (see Appendix [A]). 

In the model, we assume that the network filling factor (a„) 
is proportional to the network magnetic flux {(popen + 4'eph)- The 
proportionality coefficient il/(f>sat,n) relates the physical proper- 
ties of the network emission to the magnetic flux. However, it 
depends on the spatial resolution of the observations employed 
to compute it as well as on the estimate of the contribution of the 
bright components to the TSI. The SATIRE-T model suggests 
that the increase in the network emission since the Maunder 
Minimum is approximately 0.7 W/m ^. This value is close to 



Frohlich 



(2009 1 and em- 



the increase of 0.9 W/m^ estimated by^ 

ployed by Steinhilber et al. ( 2009] ). On the other hand, the recon 
struction of the bright component of the TSI by Wenzler et al. 
(2006) from magnetograms and continuum disk images indi- 
cates a contribution of approximately 0.2 W/m^ during the min- 
ima between cycles 21-22 and 22-23. This comparison stresses 
the need for higher resolution and more sensitive observations 
of the solar surface magnetic field, such as now being provided 
by the observations by the Helioseismic and Magnetic Imager 
(HMI) instrument on board of the Solar Dynamics Observatory 
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Fig. 4. Reconstructions of the geomagnetic field intensity. The VADM reconstructions by 



Knudsen et al. 



(|2008|l and benevey et al. 



( |2008| l are presented in panels (a) and (b). The VDM reconstructions by |Genevey et al.| ( |2008| l and |Korte & Consta ble ( 2005[ ) are 
shown in panels (c) and (d), respectively. The dashed lines represent 1-cr statistical errors. For reference, the VADM reconstruction 
by Knudsen et al. ( 2008j ) is presented in all panels. 



(SDO). Such data are expected to improve the estimate of the 
contribution of the bright components of the solar atmosphere to 
the TSI. 



3.2.2. Reconstructions for the Holocene 

Reconstructions of TSI for the whole Holocene are presented in 
Figure l6] The structure and color scheme are the same of the 
Figure p\ The first 4,500 years are covered by the reconstruc- 
tions based on the VADM values published by [Knudsen et al.| 
( |2008,) and partially by the one based on the VADM presented 
by [Genevey et al.| ( |2008| l. VADM and VDM reconstructions 
are jointly available from approximately 5000 BC to 1700 AD, 



which is the period covered by the Korte & Constable ( 2005[ l 
reconstruction. 

We evaluate the uncertainties for these reconstructions in the 
same way as we described in the previous section, by consid- 
ering three major sources. The uncertainties of each reconstruc- 
tion due to the uncertainties of the computation of the open flux 
are shown in panel (b). It is noticeable that the uncertainties 
are higher in the early part of the time series. The uncertainty 
estimated for the KN08 reconstruction decreases after approxi- 
mately 6000 BC to the level of the other reconstructions. In order 
to take into account the uncertainty related to the conceptually 
different paleomagnetic reconstructions, we consider separately 
the periods that VADM and VDM reconstructions are and are 
not jointly available. For the period that they are jointly avail- 
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Fig. 5. (a) Reconstructions of the total solar irradiance for the last 3,000 years. All reconstructions are based on the INTCAL04 
A '"*C data set, but employing different reconstructions of the geomagnetic dipole momentum. The VADM reconstructions by 
Knudsen et al . (2008) (KN 08/blue line) and [Genevey et aL] ( |2008|l (GN08-8k/thin red line ) are presented. In addition, the VDM 
reconstructions by Genevey et al. (2008 i(GN08/green line) andlKorte & Constable (2005) (KC05/gray line) are also plotted. For 



reference, 10-year running means of the SATIRE-T reconstruction are displayed from 1700 to the present (KlO/thick red line), (b) 
Estimate of the error due to the uncertainties of the open flux reconstructions. The red dotted line presents the uncertainty of the TSI 
computed using SATIRE-T. (c) Difference between the maximum and minimum values of the TSI computed employing different 
paleomagnetic models, (d) Estimated total uncertainties of the reconstructed TSI employing SATIRE-M model. 



able, i.e. after 5000 BC, we follow the method described in the 
previous section. For the period prior to jointly reconstructions, 
we linearly extrapolate the trends bound between 5000 BC and 
1700 AD. The results are presented in panel (c). The dashed 
green line in this panel displays the number of reconstructions 
employed to compute the maximum and minimum values (scale 
on the right hand side). While the uncertainty due to the pale- 
omagnetic models is on average at the same level as the uncer- 
tainty of the SATIRE-T model from 5000 BC to approximately 
1000 BC, during the last 3000 years the uncertainty is lower. The 



total uncertainty for each reconstruction, which was computed 
as presented in Sect. 3.2.1, is shown in panel (d). The overall 
uncertainty is about 0.4 W/m^. 



As can be gathered from Fig. A. 2 the TSI distributions on 
this longer time scale are significant closer to Gaussians. The 
difference to the more skewed distributions found for the last 
3000 years has to do with the irregular distributions of grand 
minima. 



12 



Vieira, Solanki, Krivova, and Usoskin: Evolution of the TSl/SSl 



- KN08 GN08-8k GNOB KC05 K10 



1366.5 




1364.5 



0.8 



2000 



0.6 



0.4 



t 0.2 



0.8 
^ 0.6 
i. 0.4 
to.2 


0.8 
^ 0.6 
i. 0.4 
to.2 




(b) 



-8000 



-6000 



-4000 
Time 



-2000 



2000 



(c) 


















i 







-8000 



(d) 



-6000 



-4000 
Time 



-2000 



2000 



Only VADM reconstructions 



-8000 



-6000 



-4000 
Time 



-2000 



2000 



Fig. 6. Reconstructions of the total solar irradiance for the Holocene. The figure has the same structure as Fig.js] but for the period 
9495 BC to the present. See the caption of Fig. |5] for more information. While the green line in panel (c) indicates the number 
of reconstructions employed to compute the uncertainty of the TSI due to the uncertainty of the evolution of the geomagnetic 
dipole moment, the black line represents the long-term trend of the uncertainty of the evolution of the geomagnetic dipole moment 
(y - -1.7 ■ 10^5 * f -f 0.15 [W/m^]). The period for which only VADM reconstructions are available is indicated in panel (d). 



3.2.3. Comparison between TSI reconstructions based on 
'"^C and "'Be records. 

In this section, we comp are the TSI reconstruct ion based on 
'"Be records presented by Steinhilber et al. ( 2009 i and two TSI 
reconstructions based on C computed in this paper. Figure 
shows the TSI reconstructions by Steinhilber et al. ( 2009| l 



(S09; VADM; green Une) and the reconstructions based on the 
Knudsen et ai:](|2008l l (KN08; VADM; blue line) and |Korte &| 



Constable| ( |2005| l (KC05; VDM; gray line) paleomagnetic mod- 
els since approximately 9500 BC. For reference, 10-year run- 
ning means of the SATIRE-T reconstruction are displayed from 
approximately 1640 AD to the present (KIO, thick red line). The 
S09 reconstruction, which was obtained from 40-year running 



means of "'Be, was resampled to decadal values employing a 
linear interplolation. An offset of 1365.6 W/m^ was introduced 
to match the present value of the TSI (PMOD composite). 

The TSI reconstructions based on different approaches are 
in the same variability range, while producing some discrep- 
ancies at mid to long-term time scales. Although the values of 
the three time series lie in the same interval and display simi- 
lar trends prior to 1000 BC, the reconstructions based on '^C 
and '"Be diverge during the last 3000 years. Figure |7j3 offers a 
better visualization of the structural difference between the time 
series. In this figure, the differences between 120-year running 
means of the models based on '"'C and '"Be records are pre- 
sented. For the period analyzed, the standard deviation of the 
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differences KN08-S09 and KC05-S09 are approximately 0.32 
and 0.34 W/m^, respectively. Note that the differences between 
'"^C and ''^Be reconstructions are larger than the differences due 
to the uncertainty of the geomagnetic field. Panels (c) and (d) 
of Fig. |7]show the scatter plots of the S09 reconstruction versus 
the KN08 and KC05 reconstructions, respectively. The correla- 
tion coefficients are indicated in the panels. As a consequence 
of the structural differences of the time series, the correlation 
coefficients between the reconstructions based on '"^C and '"Be 
records are relatively low. 

During the last 3000 years the S09 reconstruction suggests a 
higher level of TSI compared to the KN08 and KC09 reconstruc- 
tions. For a better visualization, we present in Fig.|8]the evolu- 
tion of the time series from 500 AD to the present. Near the end 
of the Maunder Minimum, the difference is approximat ely 0.5 
W/m^ if we require the TSI reconstruction by Steinhilber et al 



( 2009 | l to match the present value of the TSI. Since the two re- 



constructions are based on different approaches and cosmogenic 
isotopes, it is difficult to identify the major source of the dis- 
crepancy. Different long-term trends in the and '°Be data 
certainly play an important role ( jVonmoos et al.|2006[ |Usoskin| 
al.|2009 1, but methodological differences may also contribute. 
Note that the TSI reconstructed by S09 also shows a much flatter 
trend than that obtained by Krivova et al.| ( 2010a, ) from telescop- 
ically measured group sunspot number (KIO). 



3.2.4. Discussion of examples of reconstructed irradiance 
trends 

The reconstructions indicate that approximately 2450 years ago 
the TSI was at the same level as is observed at the present time. 
Figure |9] shows the TSI reconstructions from 1100 BC to 100 
AD. According to the VDM and VADM reconstructions, the 
peak value ranges from approximately 1366.0 to 1366.2 W/m^, 
respectively. This peak was reached about 330 years after the 
preceding grand minimum, which was observed in 765 BC. 
The TSI values during this grand minima computed employing 
different paleomagnetic reconstructions ranges from 1364.6 to 
1364.8 W/m^. This grand minimum was previously discussed 
by [Eddy] ((T977bla| i, |Stuiver|([T980)l, IStuiver & BraziunasI i 



|Goslar| ( |2003| ) and lUsoskin et al.| ( |2007| l. For reference, the 



1989) 1, 



decadally averaged SATIRE-T reconstruction (thick red line) is 
shifted 2455 years to match the minimum value observed in 765 
BC. The dashed blue line connects the beginning and the end 
of the SATIRE-T reconstruction. The amplitude and time scale 
of the TSI evolution from the previous grand minimum, which 
occurred at approximately 765 BC, to the grand maximum is 
similar to the amplitude and time scale reconstructed from the 
end of Maunder Minimum to the present. However, there are 
also noticeable differences of the evolution of the two periods. 
For roughly 200 years after the 765 BC grand minimum the TSI 
is higher than for the same interval after the Maunder Minimum. 
This maximum around 430 BC was followed by the next grand 
minimum, which is seen as a large and very rapid decrease of 
the TSI. This extended minimum lasted for approximately 150 
years. The minimum value of the TSI, which ranges from 1364.8 
to 1365.0 W/m^, is reached in 345 BC. Note that afl 3 plotted 
grand minima have the same TSI level within the uncertainty 
limits. The similarities between the evolution of the TSI for the 
ancient and the present periods cannot be used to predict the oc- 
currence of a grand minimum in the near future, but shows that 
very rapid and deep drops in the irradiance can occur starting 
from a high level of TSI, such as the present. 



The reconstructions also suggest that from about 3000 BC 
to 2000 BC the peaks in solar irradiance exhibited values at a 
similar level as observed at the present time. Figure 10 shows 
the reconstructions and uncertainties from 3000 BC to 1000 BC. 
As discussed in the previous subsection, the TSI reconstruction 
based on the VDM reconstruction ( |Korte & Constable|2005| ) dis- 
plays higher values, which are for some periods slightly higher 
than the averaged value observed during the last three solar cy- 
cles. The structure of the signal during this period is clearly 
different from the episode described in the previous paragraph. 
With the exception of a period close to 2860 BC, pronounced 
grand minima are not observed, although large amplitude oscil- 
lations are present. Obviously, the TSI exhibited a rather differ- 
ent behavior during that period of time than in the last 400 years. 

3.3. Construction of a TSI time series for investigations of 
the climate evolution during the Holocene 

One of the main motivations for reconstructing the TSI over 
the Holocene is to produce a consistent record of the exter- 
nal heat input into the coupled atmosphere/ocean/land system 
that is needed to investigate the evolution of the Earth's climate 
employing general circulation models ( Schmidt||2010[ [Schmidt| 
et al. 2010). However, in the previous sections we have produced 
multiple TSI time series whose validity and length depends on 
that of the particular geomagnetic field reconstruction underly- 
ing it. For use with climate models, however, it is desirable to 
construct an optimal time series extending through the whole 
Holocene. In this subsection, we describe a procedure for com- 
bining the SATIRE reconstructions during different periods of 
time. This is not completely trivial due to the different lengths 
of the time series we are averaging and the fact that they partly 
have different median values. 

In order to produce a representative time series, we divide 
the Holocene into four periods according to the available data 
for each period. Table |4] summarizes the periods considered and 
the data available for each period. After the end of the Maunder 
Minimum we compute daily values, while prior to the end of 
the Maunder Minimum we compute 10-year averages. For the 
period for which both solar disk magnetograms and contin- 
uum images are available (period 1) we employ the SATIRE-S 
reconstruction. SATIRE-T reconstruction is used from the be- 
ginning of the Maunder Minimum (-1640 AD) to 1978 AD. 
Prior to 1640 AD reconstructions based on cosmogenic isotopes 
are used. As VDM paleomagnetic reconstructions are available 
since approximately 5000 BC, we have to treat separately pe- 
riods 3 and 4. We emphasize that the reconstructions based on 
different proxies have different time resolutions. 

For period 3, VADM and VDM reconstructions are avail- 
able. The TSI reconstructions based on conceptually different 
paleomagnetic models differ significantly from each other As 
discussed in Sect 3.2, the paleomagnetic reconstructions based 
on the VADM and VDM models provide the upper and lower 
limits of the evolution of the geomagnetic field, respectively. 
Consequently, we need to combine them to produce a represen- 
tative reconstruction for this period. As we have discussed previ- 
ously, the median values of the reconstructions for the last 3000 
years that belong to a given group are not significantly different. 
In this way, assuming that this result holds for the whole period 
starting in 5000 BC, we can compute the final TSI at time step j 
as an average of the two methods VADM and VDM taking into 
account that they represent lower (VDM) and upper (VADM) 
bounds to the true dipole moment. We do not claim that the av- 
erage of the two groups is the true value, but since the true value 
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Fig. 7. (a) Comparison between the TSI reconstruction by 

constructions obtained in this paper based on VADM and VDM paleomagnetic reconstructions by Knudsen et al.|(|2008P (KN08; 
VADM; blue line) and Korte & Constable (2005) (KC05; VDM; gray line), respectively. For reference, lO-year running means of 
the SATIRE-T reconstruction are displayed from aprroximately 1640 to the present (KIO, thick red line), (b) Differences between 
120-year running means of the models based on '^C and "'Be records, (c) Scatter plot TSI-S09 versus TSI-KN08. (d) Scatter plot 
TSI-S09 versus TSI-KC05. 



lies between these bounds, if (the true value) is more likely to lie 
close to the average value than to either of the bounds. Here we 
employ the reconstructions based on the KC05 and KN08 dipole 
moments that are representative of the VDM and VADM groups, 
respectively. 

For period 4, only VADM reconstructions are available. Here 
we employ the KN08 reconstruction as the reference. To be con- 
sistent with other periods, it is necessary to introduce an offset 
that can be estimated from the difference of the mean values of 
the KN08-VADM and KC05-VDM reconstructions for the pe- 
riod 3. In this way, we implicitly assume that the VDM paleo- 
magnetic models would have lower values than the VADM pa- 



leomagnetic models in period 4 as they have had in more recent 
times as well. The numerical value of the offset is approximately 
0.05 W/m^. In this way, discontinuities are not introduced in the 
interface between periods 3 and 4. Note that the upper bound 
limit in Fig. 11 is obtained by adding 2 ■ offset to the VADM re- 
construction. For more recent times the SATIRE-T and S ATIRE- 
S based reconstructions are appended to the '"^C based TSI re- 
constructions (see Table Hh. The final TSI reconstruction is pre- 



sented in Figure 1 1 The periods considered and listed in Table 
|4] are indicated in the figure. To avoid overloading the figure, 
10-year averages of the SATIRE-T model are plotted. 
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Fig. 8. Comparision between the TSI reconstructions based on '"^C and '"Be from 500 AD to the present. Except for the considered 
time interval the structure of the figure is identical to panel (a) of Fig. |7] 
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Fig. 9. Comparison between the evolutions of the TSI after the grand minimum observed in approximately 765 BC and after the 
Maunder Minimum. The figure follows the format of Fig. [5^. The 10-year averaged SATIRE-T reconstruction (thick red line) 
is shifted in time by 2465 years in order to match the minimum value observed in 765 BC. The dashed blue line connects the 
beginning and the end of the SATIRE-T reconstruction and emphasizes the roughly linear rise in irradiance sind the end of the 
Maunder Minimum. 



4. Concluding remarks 

In order to understand the evolution of the Earth's climate during 
the Holocene the variability of the internal and external drivers 
needs to be quantified. Among the external drivers, the solar ir- 
radiance plays a prominent role. Here, we have reconstructed 
the TSI over the Holocene. This is an extension of the SATIRE 
family of models to a longer time scale. It represents the first re- 
construction of TSI during pre-telescopic times that is based on 
a physically consistent treatment, which builds on the assump- 
tion that the magnetic field is responsible for all TSI changes. 
This assumption implies that our model cannot handle changes 



in TSI that do not originate in surface magnetism (e.g. convec- 
tive efficiency changes unrelated to surface magnetism or effects 
associated with r modes, i.e. oscillatory modes related with ro- 
tation, or changes associated with magnetic fields buried deep 
in the convection zone). On the other hand, it allows us to rigor- 
ously derive the proper relationship between irradiance and mag- 
netic flux. Earlier reconstructions of TSI based on multi-decadal 
averages of cosmogenic data (e.g. Steinhilber et al. 2009 | l are 
based on linear extrapolations of an empirical relationship be- 
tween open flux and TSI. However, our computations do not 
fully support a simple one-to-one relationship between TSI and 
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Fig. 10. Evolution of the TSI during the period of high solar activity from 3000 BC to 1000 BC. The figure follows the format of 
Fig.|5}i. 

Table 4. SATIRE Models employed to construct a representative TSI time series for the Holocene. 



Index 


Period 




Cadence 


Model 


Input Data 


1 
2 
3 
4 


1974 AD - 
1640 AD- 
5000 BC - 
9495 BC - 


2010 
1974 AD 
1640 AD 
5000 BC 


Daily 
Daily 
Decadal 
Decadal 


SATIRE-S 
SATIRE-T 
SATIRE-M 
SATIRE-M 


Disk magnetograms and continuum images 
Sunspot observations 

Cosmogenic isotopes ( '''C; VADM and VDM geomagnetic field) 
Cosmogenic isotopes ( '''C; VADM geomagnetic field) 



solar open magnetic flux on decadal time-scale, since we find 
that magnetic flux from two cycles contributes to the TSI from 
one cycle. The reconstructions carried out here are based on this 
consistently derived relationship. Another major difference be- 
tween the present and previous reconstructions of TSI over the 
Holocene lies in our careful assessment of the uncertainties in re- 
constructions of the geomagnetic field by considering the avail- 
able geomagnetic field reconstructions that reaches back at least 
3000 years. The various reconstructions of the geomagnetic field 
can lead to significant differences in TSI at earlier times. 

The measured solar irradiance can be reproduced on the ba- 
sis of the assumption that its variability is related to the evolu- 
tion of the magnetic features lying on the solar surface ( Krivova| 
et al.'2003; Wenzler et al.|20 0?). However, the spatial and tem- 
poral resolution of records of magnetic features on the solar disk 
degrades progressively as we go back in time. For example, the 
solar irradiance can be reproduced with high accuracy from mag- 
netograms and continuum images of the solar disk employing a 
physics-based model (SATIRE-S |Krivova et al. 2003[ l. However, 
these records are jointly available only for the last three solar cy- 
cles. To deal with this limitation, models of the evolution of the 
solar magnetic field can be employed to replace the direct ob- 
servations of the evolution of the magnetic features on the solar 
disk. Most of these models are based on the records of sunspot 
number, which are available since approximately 1610 AD (e.g., 
SATIRE-T (Balmaceda 2007; Ki-ivo va et al.||2007| |20I0a|l and 



other models (Lean et al. 1995; Solanki & Fligge|| 



999[ IWangI 



[eral.,2005; Tapping et al. 2007 , Crouch et al.||20d8| l), but con 
tain no information on the spatial distribution of the magnetic 
features on the solar surface and suffer from loss of temporal 



resolution at earlier times. On longer time-scales, the evolution 
of the solar magnetic flux can be evaluated based on records of 
cosmogenic isotopes stored in natural archives, although gener- 
ally at the cost of a further loss of temporal resolution. Here, we 
have shown that 10-year averages of the total and spectral irradi- 
ance can be obtained by employing the magnetic flux computed 
from cosmogenic isotopes as the input of the SATIRE model 
(SATIRE-M, where 'M' stands for millennial time scales). Note 
that the core of the irradiance reconstruction method is the same 
for all variations of the SATIRE model. The solar atmosphere is 
divided into components: sunspots, umbrae and penumbrae, fac- 
ulae, network and the quiet Sun. The evolution of the spectrum 
is described by the filling factors on diff'erent parts of the solar 
surface of each component and its time independent brightness. 
Only the evolution of the filling factors of the surface compo- 
nents is computed according to the available data. 

The reconstructed 10-year averaged TSI displays variations 
over a range of approximately 1.5 W/m^. This is slightly larger 
than the 1.3 W/m- between the Maunder Minimum and the cur- 



rent grand maximum found by Krivova et al. ( 2010a i, mainly 



because an earlier Grand minimum (around 7500 BC) displays a 
lower TSI than the Maunder Minimum and the Grand Maximum 
at the beginning of the time series displays a larger TSI than 
the just ending Grand Maximum. We point out that this range 
is uncertain as there are significant uncertainties in the compu- 
tation of the open flux at the earlier times. The major sources 
of uncertainty of the solar irradiance reconstructions during the 
Holocene come from the simplifications employed to derive the 
evolution of the magnetic flux emerging in AR and ER based 
on decadal values of the open flux and from the evolution of 
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Fig. 11. TSI weighted reconstruction since approximately 9500 BC. In order to provide a better visualization, the evolution since 
1000 BC is displayed in panel (b). The filled gray band represents region limited by the KN08-VADM and KC05-VDM reconstruc- 
tions. For reference, the red lines represent the 10-year averaged reconstruction by Krivova et al. ( 2010a| l. 



the Earth's magnetic field. In order to evaluate the uncertain- 
ties due to the evolution of the Earth's magnetic field, we com- 
pared an ensemble of reconstructions that are based on two con- 
ceptually different approaches to computing the evolution of the 
Earth's magnetic dipole. For the period that all reconstructions 
are jointly available (i.e. during the last 3,000 years) the two 
groups show significantly different median values, reflecting the 
fact that reconstructions of the virtual dipole moment may under- 
estimate the true dipole and, consequently, overestimate the solar 
parameters. The opposite is expected for the virtual axial dipole 
moment reconstructions. In addition to these uncertainties, we 
recall that the reconstructions are based on the assumption that 
carbon cycle parameters have not changed significantly through- 
out the Holocene, which may not be exactly vahd, particularly 
in the early Holocene, as discussed previously by |Solanki et al.| 



( |2004l i, lElhson et al.| ( |2006| , [Usoskin et al^ ( |2007l |2009] l, and 
others. We stress that the TSI reconstructions presented in this 
paper a re bas ed on one individual record of '"^C, and, as shown 



in Sect. 3.2.3| using Be records would probably lead to another 



long-term behaviour 

We have also produced a reconstruction that is a careful 
combination of the various TSI reconstructions we have carried 
out for the period prior to the Maunder Minimum, i.e. before 
1640 AD. Between 1640 AD a nd 1974 AD it correspo nds to 
the S ATIRE-T reconstruction of [Krivova et al.| ( [2010a| l, while 
after that we employ a SATIRE-S reconstruction based on Kitt 
Peak (Wenzl er et~aL] [2006) and MDI (following 'Ki'ivov a et aL] 
2003) magnetograms. This combined reconstruction, which is 
our best estimate of the TSI during the Holocene, is avail- 
able as electronic material to this paper and can also be down- 
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loaded from the website: | http://www.mps.mpg.de/projects/sun-| 
[clim ate/data.html 

The present analysis is in agreement with earlier reconstruc- 
tions of the solar activity during the Holocene regarding the oc- 
currence of extended periods of low and high solar activity (e.g. 



the alternative that they do not have equal medians. Note that 
the distributions need not necessarily be normal. Table A.l sum- 



Usoskin et al. 2007 1. These coincide with periods of low and 



high TSI. We noted some similarities between the evolution of 
the solar activity from the Maunder Minimum to the present 
and an episode observed approximately 2800 years ago. Both 
episodes displayed changes of the solar irradiance of about 1.3 
W/m^ in an interval of approximately 300 years. The ancient 
episode ended rather abruptly in a grand minimum, which em- 
phasizes that a very rapid drop in the irradiance can occur start- 
ing from a high level of TSI, such as the present one. However, 
we stress that we cannot employ these similarities to predict the 
occurrence of a grand minimum in the near future. 

Appendix A: Statistical properties of the TSI 
reconstructions 

Panels (a) to (d) of Fig. [AT| display the probability distributions 
of the 10-year averaged TSI values. Each panel shows the prob- 
ability distribution functions (PDF) for the normal distribution 
computed employing the respective mean and standard devia- 
tion values that are also given in each panel. In order to verify 
if the reconstructed TSI values are normally distributed, we per- 
formed a test of the default null hypothesis that the values come 
from a normal distribution, against the hypothesis that they do 
not come from this distribution. The test statistic is 



T — max \abs {ecdfiTS I) — ncdf(TS I))] 



(A.l) 



where ecdf(TS I) is the empirical cumulative distribution func- 
tion (CDF) estimated from the TSI distribution and ncdf(TSI) 
is the normal CDF with mean and standard deviation equal to the 
mean and standard deviation of the TSI distribution. Recall that 
the value of the ecdf at a given value of the TSI is obtained from 
the PDF by integrating the latter from its lower end up to that 
particular TSI value. The T values are indicated in the respec- 
tive panels. The null hypothesis is rejected at a given significance 
level when T is greater than a critical value. The test employed 
here uses a table of critical values computed using Monte-Carlo 
simulations, where the critical value for the significance level 
of 0.05 is 0.0547. We found that the null hypothesis of all re- 
constructions is rejected at this significance level, i.e., the distri- 
butions do not correspond to normal distributions. This is sup- 
ported by a visual inspection which reveals the strong asymme- 
try in all the the empirical distributions. Besides the main peak 
located above 1365 W/nr there is also a lower (and statistically 
less significant) secondary peak below 1365 W/m^, which cor- 
responds to the relatively common grand minima in this period 
of time. 

Figure |A.l^ compares the median values (red lines). 
Additionally, the upper and lower boundaries of the boxes dis- 
play the respective quartile values. The dashed lines show the 
regions extending from the boundaries of each box to the most 
extreme values within 1.5 times the interquartile range. As ex- 
pected, the VDM reconstructions display higher median values 
than the VADM reconstructions. Furthermore, as also suggested 
by the empirical CDFs, the VDM distributions are shifted in re- 
lation to the VADM distributions, and also display a larger range 
of irradiance values. 

We employ the Wilcoxon test to verify if two reconstructions 
come from identical distributions with equal medians, against 



marizes the results. The reconstructions from the same group 
(VADM or VDM) have the same median values. However, the 
reconstructions belonging to different groups disagree from each 
other. 

The statistical properties of the TSI reconstructions for the 
whole Holocene based on the paleomagnetic reconstructions 
by [Knudsen et al.| ( [SOOS] ), [Genevey eTaL| ( |2008] l, and |Korte &" 



Constable ( 2005 1 are presented in panels (a)-(c) of Fig. A. 2 
spectively. As in Fig 



A.l 



the red lines represent the probability 
distribution functions for the normal distribution. Note that each 
time series covers a different time span (see Table [3|. Tests of 
the null hypothesis that the values come from a normal distribu- 
tion, as discussed in Sect. 3.2.1, indicate that it is not possible to 
reject this hypothesis at the significance level of 0.05 for the re- 
constructions based on the Gene vey et al.| (|2008 ) paleomagnetic 
reconstruction. A visual inspection shows that the reconstruction 
based on the paleomagnetic model by Korte & Constable ( 2005 | l, 
that extends from 5000 BC to 1700, gives a slightly asymmetri- 
cal distribution. This behavior reflects the more frequent occur- 
rence of grand minima in the later part of the time series, as can 
be gathered by comparing with the far more asymmetric distri- 
butions in Fig. A. 1 
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KN08 


GN08-8k 


GN08 


KC05 


KN08 


X 




rej. 


rej. 


GN08-8k 




X 


rej. 


rej. 


GN08 


rej. 


rej. 


X 




KC05 


rej. 


rej. 




X 



KN08 



GN08-8k 




1364 1364.5 1365 1365.5 1366 1366.5 
TSI [W/m^] 
GN08 




1364 1364.5 1365 1365.5 1366 1366.5 
TSI [W/m^] 



10 



Mean : 
Median : 
Std : 



1365.30 
1365.35 
0.31 
0.000 





(b) 











1364 1364.5 1365 1365.5 1366 1366.5 
TSI [W/m^] 



KC05 



40 

B ^° 

0) 

i3 20 

O 

d 

^ 10 



Mean : 1365.41 
Median: 1355.45 
Std : 0.34 




1364 1364.5 1365 1365.5 1366 1366.5 
TSI [W/m^] 



1366 



1365.5 



1365 



1364.5 



_(e) 



L J 



KN08 



GN08-8k GN08 
Model 



KC05 



Fig. A.l. Comparison between the statistical properties of the TSI reconstructions for the period 1000 BC to 1700 AD. Panels (a) to 
(d) present the distributions of the TSI reconstructions. The mean, median and standard deviation values are indicated in the panels. 
The red lines represent the probability distributions functions for the normal distribution. Panel (e) allows a visual comparison 
between the distributions. The median (red line), lower and upper quartile values of the distributions (lower and upper ends of the 
blue boxes) are plotted in each box. The black lines extend from the end of each box to the most extreme values within 1 .5 times 
the interquartile range. Note that the interquartile range, which is also a measure of the statistical dispersion of the data, is equal to 
the difference between the third and first quartiles. 
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